﻿ Detectarea schimbării dispersiei Distincția între două ipoteze poate fi dată de schimbarea dispersiei unei statistici gaussiene Considerăm că datele gaussiene x[n] sunt statistic independenteși identic distribuite (IID) Repartițiile datelor, în cele două ipoteze, sunt 2 =∼HN () 00: [ ] 0, ; 0, 1, , N-1xn nσ 222 H σσ=>N ∼ () 11 10: xn nσ [ ] 0, ; 0, 1, , N-1 Testul NP decide că ipoteza adevărată este H1, dacă raportul de plauzibilitate depășește pragul testului, γ 1N− 211⎧⎫ () ⎨⎬ ∑ /22exp [ ]Nxn− nσ=02122; ⎩⎭ ()xH () 11pπσ == > Lγ x () 1N− () 2;11p⎧⎫ 0xH () ⎨⎬ ∑ /22exp [ ]Nxn− 02022nσ= ⎩⎭ () 0πσ Efectuăm calculele și obținem N 122N− ⎛⎞ ⎛⎞ 211 1σ ⎪⎪ 0⎧⎫ −> () ⎨⎬ ∑ ⎜⎟ ⎜⎟ 222[]exp x nγ σ 01012nσσ= ⎪⎪ ⎝⎠ ⎝⎠ ⎩⎭ După logaritmare se determină condiția de alegere a ipotezei H1ca fiind 122N− σ+ () () 2102lnln1Nγσ ′ >=> () ∑ 22[] 0xnγ ()() 011nNσσ= 01− 1 Statistica testului pentru detectarea schimbării dispersiei este deci estimatorul acesteia 1N− 21 x []Txn= ()() ∑ 0nN= În mod evident, pentru statistica testului se poate lua estimatorul MVU al dispersiei 1N− 21 x []Txn= ()() ∑ − 01nN= Pentru N suficient de mare, cele două statistici sunt, practic, identice rată, dacă Pentru N=1, testul decide că H1 este ipoteza adevă 2 ′′ >> xxγγ () Cele două repartiții, corespunzătoare celor două ipoteze, precum și cele două regiuni de decizie sunt arătate în figură Regiunile de decizie au formele ′′ : xγγ , , ∈−∞ − ∞ ≡ ∪H ()() 11Decidedacă ′′ : xγγ , ∈− ≡H () 00Decidedacă *** Acest comportament este determinat de existența unei statistici suficiente pentru parametrul în cauză Remarcăm că, atât în cazul detectării unui nivel continuu (din zgomot) cât și în cazul detectării schimbării dispersiei, trebuie să decidem între două ipoteze ale căror densități de probabilitate au anumiți parametri de valori diferite, în cele două ipoteze Pentru luarea deciziei, se estimează parametrul care se modifică și se compară valoarea estimată cu un prag În cazul general, se observă datele x care au o densitate de probabilitate parametrizată de θ, p(x;θ) Trebuie să efectuăm un test pentru parametrul θ, test descris prin θ=H 00: θ θ= H 11: θ Dacă există o statistică suficientă pentru parametrul θ, T(x), conform teoremei de factorizare Neyman-Fisher, densitatea de probabilitate a datelor se poate pune sub forma ; ( ); ( )pgT hθθ=xxx ()() Aplicând testul NP, vom decide că H1 este ipoteza adevărată, dacă xx ()() pgTθθ ( );11; Lγ == >x () xx ()() pgTθθ ( );00; Menționăm că factorul h(x) s-a simpificat Datele x sunt aceleași dar nu știm din care repartiție provin Inlocuim aceleși date în ambele modele de repartiție și vedem în care ele “se potrivesc mai bine” 2 Dacă parametrul θ este o medie, atunci aplicând teorema de factorizare avem 1N− 211⎧⎫ x =−− ()() ⎨⎬ ∑ exp [ ]Npxnθθ/22; ⎩⎭ 0222nσ= πσ () Tx () ⎧⎫ ⎛⎞ ⎪⎪ ⎜⎟ 112NN−− ⎧⎫ 211 1NNθθ ⎪⎪ =− − () ⎨⎬⎨ ⎬ ∑∑ /2222 exp [ ] exp [ ]Nxn xn⎜⎟ σ ⎜⎟ ⎩⎭ 002222nnNσσ== πσ ⎪⎪ () ⎜⎟ ⎪⎪ ⎩⎭ h⎝⎠x () gTθx; () ) ( Dacă parametrul θ este o dispersie, aplicând aceeași teoremă, rezultă Tx () ⎧⎫ ⎛⎞ ⎪⎪ ⎜⎟ 1N− 211N⎪⎪ ⎜⎟ =− ⋅ x ()() ⎨⎬ ∑ exp [ ] 1Npxnθ/2; ⎜⎟ 0222nhNθ= πσ x () ⎪⎪ () ⎜⎟ ⎪⎪ ⎝⎠ ⎩⎭ gTθx; () () OBSERVATIE IMPORTANTĂ x Dacă este un estimator nedeplasat al parametrului ce distinge între douăTθ () ipoteze, atunci detectorul NP va fi implementat ca un estimator al parametrului Vom aborda un exemplu în care trebuie să detectăm prezența unei componente continue, A, dintr-un zgomot, w[n], care nu este gaussian Datele, în cele două ipoteze, sunt =H 0: [ ] [ ]; 0, 1, , N-1xn wn n= =H 1: xn A wn n=+ [ ] [ ]; 0, 1, , N-1 Eșantioanele de zgomot, w[n], sunt statistic independente și identic distribuite (IID), dar nu sunt gausiene Fiind eșantioane IID putem scrie, pentru cele două ipoteze, relațiile 11NN−− xH ()() ( ) ∏∏ 0;[][]ppxnpwn== 00nn== 11NN−− =xH () ( ) () ∏∏ 1;[] []ppxn A pwn=− 00nn== Pentru precizarea ideilor, vom considera că zgomotul este o mixtură (un amestec) de de două zgomote gaussiene, unul cu dispersia 1, iar celălalt cu dispersia 2, contribuția fiecărei componente fiind de 0 5, așa cum se arată în figură 3 În cazul considerat, expresia densității de probabilitate a zgomotului non-gaussian este 22 ⎫⎧⎫ [] []11 11wn wn⎧ ()() ⎪⎪⎪⎪ [] exp exppwn=⋅ − +⋅ − () ⎨⎬⎨⎬ 222424ππ ⎪⎪⎪⎪ ⎩⎭ ⎩⎭ () 2[][]pwnpwn () 1 rată, dacă raportul de plauzibilitate Detctorul NP va decide că H1 este ipoteza adevă depășește pragul γ 1N− []pxn A− () ∏ xH () 10;np= ==> x () 1NLγ − () 0;pxH []pxn () ∏ 0n= Substituind expresia densității de probabilitate în raportul de plauzibilitate, obținem ⎡⎤ 22⎢⎥ 1N− ⎧⎫ −− [] []11 11wn A wn A⎧⎫ ()() ⎪⎪ ⎪⎪ ⎢⎥ +⋅− exp exp⋅− ⎨⎬ ⎨⎬ ∏ ⎢⎥ 0222424nππ = ⎪⎪ ⎪⎪ ⎩⎭ ⎩⎭ ⎢⎥ () 2[][]pwnpwn⎢⎥ () 1 ⎣⎦ γ> ⎡ ⎤ ⎥ 22⎢ 1N− ⎧⎫ [] []11 11wn wn⎧⎫ ()() ⎪⎪ ⎪⎪ ⎢ ⎥ +⋅− exp exp⋅− ⎨⎬ ⎨⎬ ∏ ⎢ ⎥ 222424ππ 0n= ⎪⎪ ⎪⎪ ⎩⎭ ⎩⎭ ⎢⎥ () 2[][]pwnpwn⎢⎥ () 1 ⎣ ⎦ Analiza acestei ultime relații pentru raportul de plauzibilitate ne arată că nu putem face vre-o simplificare notabilă, și nici nu observăm existența unei statistici suficiente Problema nu poate fi deci “redusă”; trebuie să substituim valorile măsurate pentru datele x în expresie, să efectuăm calculele și să comparăm rezultatul cu valoarea praguluiγ *** Caracteristicile de operare ale unui receptor (ROC) Performanțele unui detector de tip Neyman-Pearson (NP) pot fi redate prin grafice care dau dependența dintre probabilitatea deciziei corecte (puterea testului) și probabilitatea alarmei false (semnificația testului) Aceasta din urmă, ca valoare impusă, se reprezintă pe axa orizontală Reluând rezultate anterioare, pentru detecția unei componente continue, afectată de un zgomot alb, gaussian, dacă avem o secvență de N date, avem ⎛⎞ ′ γ == PQα ⎜⎟ 2FA⎜⎟ Nσ ⎝⎠ 2 ⎛⎞ ′ − − 122ANAγ −= ⎜⎟ () ) ( DFAPQ QQP d d== 2; 2 ⎜⎟ σ Nσ ⎝⎠ Relația se reprezintă grafic în figura următoare, pentru o valoare dată a raportului semnal/zgomot, SNR (sau ENR, dacă preferați) și constituie ROC (Receiver Operating Characteristic) 4 Dacă punem pragul la minus infinit, se detectează toate țintele, dar și toate alarele false →−∞ = = Dacă , atunci 1 si 1 (punctul A pe figură)PPγ DFA′ Punând pragul la plus infinit, na apare nici o alarmă falsă, dar nici nu se detectează ținte →+∞ = = Dacă , atunci 0 si 0 (punctul B pe figură)PPγ DFA′ Curba este trasată pentru d=1 Dacă probadilitatea alarmei false este 0 4, rezultă o probabilitate de detecție de 0 78 Plasarea ROC deasupra bisectoarei este urmare a faptului că recurgem la un detector, adică DFAPP> Bisectoarea constituie ROC pentru un receptor care ignoră complet datele și ia decizii în mod complet aleator, “dând cu banul” Pentru un astfel de receptor am avea H === {}{} ""FAPPcap Pcap p⎫0""; ⇒ PP= DFA⎬ H = {}{} ""DP P cap P cap p==1""; ⎭ deoarece se ignoră complet condiționarea adusă de datele măsurate Se ajunge la egalitatea celor două probabilități, deci la o curbă ROC identică cu bisectoarea AB 5 În figură se arată o familie de caracteristici de operare (ROC), pentru mai multe valori ale coeficientului de deflexie, care este de fapt echivalent raportului semnal pe zgomot, SNR La aceeași valoare admisă pentru probabilitatea alarmei false, odată cu creștere valorii pentru d, crește și probabilitatea de detecție corectă Probabilitatea minimă de eroare Abordarea bayesiană În unele cazuri putem aloca probabilități apriorice ipotezelor Ne exprimăm astfel o încredere “ante-factum” în plauzibilitatea ipotezelor Acesta e, spre exemplu, cazul transmisiilor binare digitale în care cele două simboluri transmise, 0 sau1,apar cu probabilități egale cu 0 5 Celor două ipoteze 0: a fost transmis simbolul "0"H 1: a fost transmis simbolul "1"H le alocăm probabilități apriorice egale HH {}{} 0112PP== Modalitatea de abordare în care asociem ipotezelor implicate în problema de detecție probabilități apriorice egale, se numește abordarea bayesiană În această abordare se poate defini o probabilitate medie de erore, fără a ține seama de natura erorii ratădecide , când este adevăratăPP P=+HH HH {}{} 01 10edecide , când este adevă i atunci când ipoteza unde este o probabilitatecondiționată de a decide H {} jiPHH j Probabilitatea medie de erore se rescrie, sub forma adevărată este H HH H HH H {}{} {}{} 01 1 10 0ePP P P P=+ decide , de faptul că este adevăratăconditionatăPP=HH H H {} {} jjii ;PHH {} i când Hj este adevărată, ie între , probabilitatea de a decide H face distincțjiVom {} jiPHH fără însă a-i aloca oprobabilitate apriorică, și probabilitatea condiționată Un detector care decide că ipoteza H1dacă H {} {} 10pPxH : )(1γ>= H decide 1Se dacă {} {} 01pPxHH 6 minimizează probabilitate medie de eroare, Pe Raportul de plauzibilitate condiționată, se compară cu un prag γ fix, ce rezultă din raportul probabilităților apriorice ale ipotezelor, luate în succesiune inversată Se practică și notația 1H H {} {} 10pPxH γ=> {} {} 01pPxHH sau una mai pretențioasă, ce indică modul în care se decid ambele ipoteze Dacă, așa ca în transmisiile binare, cele două probabilitați apriorice sunt egale, adică HH {}{} 0112PP== pragul testului bayesian de ipoteze este egal cu unu și deci decizia se ia conform cu 1H =>xxHH {}{} 10; 1 (2)ppγ Datele x fiind recepționate, deci fixate, expresiile {}{} ppxxHH10 si sunt, de fapt, plauzibilități condiționate Din regula de decizie se vede că decizia se ia în favoarea ipotezei cu probabilitatea condiționată maximă Detectorul bayesian, pentru cazul probabilităților apriorice egale, se numește detector de plauzibilitate (condiționată) maximă Detectorul se notează cu acronimul ML (Maximum Likelihood Detector) Vom exemplifica cu un detector de probabilitate minimă de eroare, pentru detecția unui nivel continuu, A>0, dintrun zgomot alb, gaussian Problema de detecție este =H 0: [ ] [ ]; 0, 1, , N-1xn wn n= =H 1: [ ] [ ]; 0, 1, , N-1; A>0xn A wn n=+ 2 , σw0I∼ () uN O problemă de acest fel apare în transmisiile OOK, când semnalele utilizate sunt = 01[ ] 0; [ ]snsnA= Cum probabilitățile apriorice sunt egale, receptorul (detectorul) bayesian ce minimizează probabilitate (medie) de eroare decide conform regulii 1N− 211⎧⎫ exp [ ]xn A−− () ∑ H 1/222N⎨⎬ 022nσ= ⎩⎭ πσ () => 11Nγ− 211⎧⎫ () ⎨⎬ ∑ /22exp [ ]Nxn− 0222nσ= ⎩⎭ πσ () După simplificare și logaritmare rezultă 1N− 1⎛⎞ +> ∑ ⎜⎟ 22[] 0AxnNA−− 02nσ= ⎝⎠ În final, detectorul bayesian ce minimizează, în acest caz, probabilitatea de eroare este 1H 1N− 1A > []xxn= ∑ 02nN= Pentru a determina valoarea erorii medii minime, avem în vedere că media eșantion este repartizată, în cele două ipoteze, conform cu 7 2 ⎧ σHN () Nconditionată de ipoteza 00, ⎪ x⎨ ∼ 2 σH N () ⎪ AconditNionată de ipoteza 1, ⎩ Se aplică regulile de calcul stabilite până acum și obținem +⋅HH HH {}{} 01 10e12 12PPP=⋅ ⎡⎤ ⎢⎥ ⎧⎫ ⎧⎫ 1⎪⎪ ⎪⎪ ⎢⎥ HH +> ⋅ ⎨⎬⎨⎬ 2 2Px A Px A= ⋅HH {}{} 1 2 2Px A Px A=−10 ⎜⎟ 2⎢⎥ ⎢⎥ ⎝⎠ ⎣⎦ ⎛⎞⎛⎞ 220AA A⎡⎤ −− ()() 1 1QQ⎢⎥ ⋅ =− + ⎜⎟⎜⎟ 22⎜⎟ ⎜⎟ 2 ⎢⎥ NNσσ ⎝⎠⎝⎠ ⎣⎦ ⎧⎫ ⎛⎞⎛⎞ 221AA⎡⎤ ⎪⎪ 1 1QQ⎢⎥ =−− + ⋅ ⎜⎟⎜⎟ 22⎨⎬ ⎜⎟⎜ ⎟ 2NNσσ ⎢⎥ ⎪⎪ ⎝⎠⎝⎠ ⎣⎦ ⎭⎩ 2 1NA⎛⎞ > ⎜⎟ 2 ; 0AQ= ⎜⎟ 4σ ⎝⎠ În concluzie, probabilitatea de eroare minimă are expresia 2 22NA ) ( PQd d== 2e4; σ și nu mai depinde de nimic altceva, decât de coeficientul de deflexie (adică de SNR) *** O altă formă a detectorului bayesian ce minimizează probabilitatea (medie) de eroare se poate deduce din regula (1), rescrisă sub forma 1H xxHH HH {}{} {}{} 11 0 0 ( 3pP p P>) Reamintim formula lui Bayes {} {} iipPxHH ; 0, 1Pi== {} ixH px () care leagă probabilitatea datelor condiționate de ipoteză, P{xIH}, de probabilitatea ipotezei condiționată de date, P{HIx}, numită și probabilitate aposteriori În formula lui Bayes cu p(x) s-a notat expresia xx xHH HH (){}{} {}{} 00 11pp P p P=+ Regula de decizie (3) devine, după aplicarea formulei lui Bayes 1H >xx xx xHH ()() () {}{} 10; 0Ppppp> Simplificând, rezultă detectorul 8 1H xxHH {}{} Pp>10 (4) Se alege deci acea ipoteză care are probabilitatea aposteriori maximă Facem câteva precizări 1) Detectorul care minimizează probabilitatea medie de eroare, Pe, pentru orice i}, se numește detector de probabilitate probabilități apriorice ale ipotezelor, P{H aposteriori maximă, MAP (Maximum Aposteriori Probability detector) și decide conform relației (4) 2) Detectorul care minimizează probabilitatea medie de eroare, Pe, atunci când probabilitățile apriorice ale ipotezelor, P{Hi}, sunt egale între ele, se numește detector de plauzibilitate maximă, ML (Maximum Likelihood detector) și decide conform relației (1) *** Pentru exemplificare, ne vom referi la cazul discutat deja, cel al detecției unei componente continue, de valoare A=2 5, dintr-un zgomot alb, gaussian, repartizat N(0, 1), detecție ce se face pe baza unui singur eșantion Dacă probabilitățile apriorice ale celor două ipoteze sunt egale, detectorul va fi unul de tip ML și se ajunge la 1H / 2 1 25xA=> Repartițiile eșantioanelor x posibileși regiunile critice sunt cele din figură Vom discuta acum acelaşi exemplu, dar pentru probabilităţi apriorice diferite, anume HH {}{} 011 5 si 4 5PP== Detectorul bayesian va fi de tip MAP, care decide că ipoteza adevărată este H1 dacă 2⎧⎫ 2 51x− () ⎪⎪ exp1− ⎨⎬ 22 pxπ H () ⎩⎭ 115⎪⎪ exp 2 5 3 125x== −= > {} 244⎧⎫ pxH () () 0 1x⎪⎪ exp5− ⎨⎬ 22π ⎪⎪ ⎩⎭ Logaritmând obţinem 3 125 ln 4− 0 696x>≅ 2 5 9 Pragul de comparare al eşantionului x a scăzut de la valoarea de 1 25, corespunzătoare detectorului ML, la valoarea de 0 696, ca urmare a probabilităţilor apriorice neegale În figură se arată repartiţiile eşantioanelor x şi noul prag (MAP) Riscul lui Bayes (sau riscul bayesian) O generalizare a abordării detecţiei, prin minimimiyarea erorii medii, ateşează fiecărei decizii ce se poate lua, un “cost”, adică o pondere a deciziei respective în bugetul erorii este costul implicat de luarea deciziei , atunci când este ipoteza adevăratăCHH ijij Dacă indicii sunt egali, decizia este una corectă, iar dacă nu sunt egali, decizia este una te un un incorectă Din considerente ce se expun în materialul scris, C10>C01 Se defineş cost mediu al tuturor deciziilor care se pot lua, ţinând seama de costurile acţiunilor întreprinse, de probabilităţile condiţionate şi de probabilităţile apriorice Avem HH H HH H {}{}{} {}{} 00 0 0 0 01 0 1 1EC C P P C P P=+ HH H HH H {}{} {}{} CP P CP P++10 1 0 0 11 1 1 1 Acest cost mediu (în sens statistic), se numeşte şi riscul lui Bayes, şi se notează cu R, şi pentru deciziile binare are expresia 11 RHHH {} {} ij i j jCP P= ∑∑ 00ij== Pentru costurile i 1CC CC==== 00 11 01 100 ş riscul lui Bayes devine chiar probabilitatea medie de eroare, prezentată pentru cazul transmisiilor binare Sunt la fel de “grave” pentru o transmisie transformarea, la recepţie, a unui unu în zero, ca şi transformarea unui zero în unu De aceea costurile deciziilor incorecte sunt egale Deciziile corecte nu dau erori şi de aceea costurile lor sunt nule În orice caz, cu în ipoteza, rezonabilă, că o decizie incorectă are un cost mai mare decât dacă nu s-ar fi produs o decizie incorectă, adică CC CC>>i 10 00 01 11 ş detectorul care minimizează riscul lui Bayes este cel ce decide că ipoteza adevărată este H1conform relaţiei 10 xH H () (){} 110 00 0pCCP− = 0 (6)γ>> xH H () {} () 01 11 10CCPp− ***Detectorul ce minimizează riscul lui Bayes (6), pentru care avem costurile i 1CC CC==== 00 11 01 100 ş 1 se face dacă alegerea ipotezei H H () {} 10pPxH 0 (7γ=>>) H {} () 10PpxH Detectorul ce minimizează riscul lui Bayes (6), pentru care avem costurileşi probabilităţile apriorice ale ipotezelor CC CC P P== == =HHi {}{} 00 11 01 10 0 10; 1 ş alegerea ipotezei H1 se face dacă () 1pxH 1γ=> (8) () 0pxH Prin urmare, detectorul (6) care minimizează riscul lui Bayes are ca şi caz particular detectorul MAP (7), care, la rândul său, are ca şi caz particular detectorul ML (8) Riscul lui Bayes sau costul mediu statistic al tuturor deciziilor posibile, devine probabilitatea minimă de eroare în cazul detectoarelor MAP şi ML TESTAREA IPOTEZELOR MULTIPLE Uneori avem de ales între mai mult de două ipoteze Aşa este cazul în recunoaşterea formelor sau în telecomunicaţii, când ne servim în transmisii de mai mult de două simboluri Pentru transmitera grupurilor de câte doi biţi (dibiţi) utilizăm M=4 simboluri, iar pentru transmiterea grupurilor de câte trei biţi (tribiţi) utilizăm M=8 simboluri, etc În astfel de transmisii se apelează la minimizarea erorii medii Deşi şi abordarea Neyman-Pearson este posibilă, ea se aplică mai rar Ea este potrivită în anumite probleme de căutare în bazele de date Presupunem că trebuie să alegem între M ipoteze HH H {} 01 1, , , M− iatunci când ipoteza adevărată este Hjcostul Notând cu Cijcostul alegerii ipotezei H mediu, statistic, al tuturor deciziilor posibile este dat de riscul lui Bayes 11MM−− RHHH {} {} ij i j jCP P= ∑∑ 00ij== În cazul comunicaţiilor digitate, aşa cum am arătat, costul fiecărei acţiuni este 1 pentru ij≠ ⎧ ⎨ ijC= 0 pentru ij= ⎩ caz în care riscul lui Bayes devine chiar probabilitatea medie de eroare, Pe Urmărim să proiectăm un detector care să minimizeze riscul bayesian 11 Dintre toate cele M ipoteze o vom alege pe aceea care minimizează expresia 1M− 0, 1, , M-1 (9)CCP i==xxH () {} iijj∑ 0j= pentru toate cele M valori ale indicelui i Altfel spus, cu datele x măsurate, calculăm expresia (9) pentru cele M valori ale indicelui i, şi CCC C≤xxx xH , , , decizia este: ipoteza este adevărată ()()()() {} 01 1dacăkM− k Regula dată prin relaţia (9) şi afirmaţiile subsecvente ei, asigură minimizarea riscului lui Bayes *** Ne vom referi la comunicaţiile digitale, în care se transmit M>2 simboluri, cu costurile aferente unei transmisii digitale de simboluri aprioric echiprobabile Riscul lui Bayes =Pe Expresia (9) devine devine atunci eroarea medie, R 11MM−− 0, 1, , M-1CPPPi==−=xx xxHHH () {} {}{} ij ij∑ ∑ 00jj== ji≠ Suma marcată (cu roșu) este o constantă, odată ce datele x sunt realizate, deoarece ea nu depinde de indicele i Minimizarea expresiei Ci(x) se realizează maximizând termenul ce se scade, adică P(HiIx) Rezultă regula de decizie MAP M-ară: Dacă cea mai mare dintre toate probabilitățile aposteriori este P(HkIx), adică , se decide că este adevărată (1,0)PPik>≠xxHH H {}{} ikk 0, 1, , -1;iM= În fapt, cu datele xrecepționate, se calculează M probabilități aposteriori, anume ⎧⎫ ⎪⎪ xx xHH H {}{}{} 01 1, , , MPP P⎨⎬− ⎪⎪ ⎩⎭ Se caută cea mai mare valoare dintre aceste probabilități aposteriori și fie ea {} kPxH Atunci se ia decizia Hk , decizie care minimizează Pe, probabilitatea medie de eroare *** Ca și în cazul detecției binare, decizia se ia maximizând probabilitatea aposteriori, motiv pentru care se spune că (10) este o regulă de decizie MAP, M-ară Un caz uzual în telecomunicații este cel în care ipotezele (simbolurile) au probabilități apriorice egale, adică 1 , 0, 1, , -1PiM== {} ixH M Costurile sunt cele normale într-o transmisie digitală Aplicând formula lui Bayes avem H {} {}{} ii ipP pxxHH ; 0, 1, , M-1Pi===xH {} i pMpxx ()() =+++ HH H H xx x xHH (){} {} {} {} {} {} 1100 11 MMpp P p P p P−− Odată ce datele x sunt realizate (măsurate), valoarea lui p(x) nu depinde de indicele i, adică nu depinde de ipoteză Aceeași afirmație este valabilă pentru numitorul Mp(x) Maximul probabilității aposteriori se obține deci pentru maximul numărătorului p(xIH i) care este, de fapt, plauzibilitatea ipotezei Hi (datele x sunt fixate) 12 Rezultă, din cele prezentate, regula de decizie ML, M-ară: Dacă , 0, 1, , -1 (11ppiki M=>≠xxHH) {} {} ik se decide că Hkeste ipoteza adevărată În fapt, cu datele x recepționate se calculează cele M plauzibilități ale ipotezelor xx xHH H {}{}{} {} , , Mpp p−101, Se caută care este ceamai mare dintre cele M valori calculate, șifie ea {} kpxH Atunci se ia decizia H k, decizie care minimizează Pe, probabilitatea medie de eroare *** Ca și în cazul detecției binare, decizia se ia maximizând plauzibilitatea ipotezei, motiv pentru care se spune că (11) este o regulă de decizie ML, M-ară Observație În cazul unei detecții MAP, M-are, se urmărește maximizarea probabilității aposteriori, {} {} iipPxHH ; 0, 1, , M-1Pi== {} ixH px () Cum p(x) nu depinde de indicele i (de ipoteză), maximizarea se obține maximizând numărătorul, sau logaritmul său ln lnpp+xHH () () ii Putem enunța o formă modificată a regulii de detecție MAP, M-are: Dacă ln ln ln ln ; 0, 1, , -1 (1 )2ppppiM++=>xxHH HH () () () () iikk atunci se ia decizia Hk , decizie care minimizează Pe, probabilitatea medie de eroare *** Vom exemplifica cu un caz de detecție ML, ternară (M=3) efectuată pe baza unui singur eșantion, x , (N=1) Ipotezele diferă prin medie, zgomotul w[n] ce afectează x , este alb, gaussian Pentru N eșantioane, cele trei ipoteze sunt +H [][] 0: ; n=0, 1, , N-1; A>0xn A wn=− H [] [] 1: xn wn= ; n=0, 1, , N-1; 2 =+H ∼N [] [] [] () 2: xn A wnσ ; n=0, 1, , N-1; w n 0, Pentru detecție vom aplica regula ML (11), dar pentru un singur eșantion Ipotezele au probabilități appriorice egale, motiv pentru care putem efectua o detecție de tip ML 13 În figură se prezintă densitățile de probabilitate ale datelor, în cele trei ipoteze, penru A=2 5 și dispersia unitară Determinarea celor trei regiuni de decizie se face simplu, analizând care dintre curbele de densitate este deasupra Rezultă imediat că >⇒HH HHH [][][][][] ()()()() x R px px px px∈>00 1 0 2 00: 00 si 00decide >⇒HH HH H [] [][][][] ()()()() x R px px px px∈>11 0 1 2 10: 00 si 00decide >⇒HH HHH [] [][][][] ()()()() x R px px px px∈>22 0 2 1 20 : 0 0 si 0 0 decide Dacă dorim să aplicăm metoda ML pentru N eșantioane, N>1,se pleacă de la 1N− ⎧ 11i=20, 1, 2 ⎧⎫ − xH [] () () ⎨⎬⎨ iiNpxnA=−22exp ; 2∑ =− = = ⎩⎭ = ⎩ 0120;0;22nAAA A Aσ πσ () și se caută pentru care ipoteză apare maximul ei, dacă cele N eșantioane sunt realizate Maximizarea se produce atunci când se minimizează expresia de la exponent, adică 11NN−− 22 2 DxnA xnxxA=− −+−= () [][] ()() () ii i∑∑ 00nn== 111NNN−−− 22 = 2xnx xA xnx xA−+ − −+ − ()() [][] ()() ∑ ii∑∑ 000nnn=== 1N− 22 2xn x x A Nx Nx NxA=−+−−+− ()()() [] () ii∑ 0n= 1N− 22 xnNx Ax=−−+ () [] () i∑ 0n= Termenul marcat (cu roșu) nu depinde de ipoteză Minimul are loc pentru minimul modulului ixA− În fapt, se determină media eșantion a datelor se calculează cele trei module și se alege −≡+ − {}{} dintre , , , ,xAxAxA xAxxA−−012minimul Minimul găsit va determina și ipoteza selectată pentru minimizarea erorii medii Probabilitatea medie minimă de eroare, pentru o detecție MAP M-ară se determină cu 22 ⎛⎞ 21MdNA− ⎜⎟ PQd==e22; ⎜⎟ 4Mσ ⎝⎠ 14